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Dalitz plot of the 7/ — ► r/7r + 7r~ decay is studied using the data collected with the 
VES spectrometer in two different exclusive reactions. The coefficients for the matrix 
element squared decomposition are measured on the largest statistics of rj' decays 
reported so far. 
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1 Introduction. 



In spite of the long history, the rf — > r)7nr decay is still not well studied and 
comprehended. The interest in the given decay is stimulated primarily by the 
nature of r/-meson itself — the large gluon component of rj' -meson should have 
strong influence on the dynamics of its dominant decay rj' — > i]7i7i [1]. 

Theoretical treatment of the rj' decays is based as a rule on the effective chi- 
ral Lagrangians (ChPT). Early calculations of the rf — > r]7T + 7T~ decay matrix 
element (within the lowest order of the expansion) resulted in the width es- 
timation which was few times less than the experimental value (see ref. [2]). 
The analogous calculations in paper [3] for the Dalitz-plot parameters have 



Preprint submitted to Physics Letters B 



2 February 2008 



overestimated the slope significantly. Phenomenologicaly the mass of the rj' is 
comparable with the theory scale and this can degrade the accuracy of the 
momentum expansion. 

The rescattering in the final state could lead to significant contribution of the 
meson resonances, in particular of the scalar octet, up to their dominance ("sa- 
turation") in the coupling constants (see [4] and references therein). The high 
mass of the 7/ makes its decays sensitive to the parameters of the intermediate 
resonances in the relevant approaches. 

The Dalitz plot for decay rj' — > ??7r + 7r - is usually expressed in variables 



here T v , T n denote the kinetic energies of mesons in the 7/ rest frame, and 
Q = T v + T n + + T n - 1 . There are two different parametrizations of event 
density on Dalitz plot in literature: the PDG [5] parametrization (which is 
also called linear parametrization) 

\M\ 2 = \1 + aY\ 2 + cX + dX 2 ; (2) 

and more general decomposition 

\M\ 2 = 1 + aY + bY 2 + cX + dX 2 ; (3) 



where a is a complex parameter and a, b, c, d are real parameters. Parameter 
c = if the Charge Parity conservation holds. Parametrizations (2) and (3) are 
equivalent if b > A measurement of parameter a at the statistics of ~ 1400 
events is reported in publication [6], here the Im(a) is neglected, and Re(a) = 
—0.08 ± 0.03. Other measurement is performed by CLEO collaboration [7] at 
the statistics of ~ 6700 events and significant background, it yields Re(a) = 
-0.021 ±0.025. 

Matrix element of the 77' — >■ r/7r°7r° decay has been studied in [8] at the statistics 
of 5400 decays, the value Re(a) = —0.058 ± 0.013 has been obtained. 

Decay channel rj' — > i]7i + 7i^ has been used already in search for C-parity 
violation [9]. No statistically significant asymmetry between ir + and ir~ has 
been observed, the measured asymmetry is —0.04 ± 0.08. 



1 The reconstructed 7/ mass and the kinetic energies used, the systematics from 
experimental smearing should be corrected by Monte Carlo methods. 
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Fig. 1. Reconstructed 77-mass distributions (Real data, histograms and Monte Carlo, 
smooth curves) a) in Diffractive-like reaction; b) in charge-exchange reaction. 
Events from the rf peak are selected. Background under rf is subtracted. 

Measured form-factors of the Dalitz plot have been used in a study of hypo- 
thetical qq and qqqq nonets of scalar mesons in bag model framework [10]. 
Later on the rf decay width and form-factors have been considered in the chi- 
ral Lagrangian approach [11,12,13,14,15]. Measured values of parameter a in 
rf — ► ?77r 7r and rf — > r/7r + 7r~ decays are used in a study of possible nonet of 
scalar mesons in article [12]. Values of parameters a,b,d in the parametriza- 
tion (3) are discussed in ChPT, taking into account the final state interactions 
in paper [13]. This approach is developed in recent reports [14,15]. 

The present work extends our previous study [16]. 



2 Experimental procedure. 

This study is based on the statistics acquired by VES experiment for inter- 
actions of 7r~ beam at the momentum of 27 GeV/c on Be target, in reac- 
tions ir~p — > rfn (charge-exchange, ~ 14.6 • 10 3 decays above background 
accepted for the analysis) and tt~N — > r/'n^N (diffractive-like production, 
~ 7 • 10 3 decays). The r\ mesons off rf decays were detected in rj — > 77 mode 
2 . The background level within the mass band of ±21 MeV near the rf mass 
was estimated as ~ 5% for the charge-exchange channel and ~ 15% for the 
diffractive-like channel. Standard background subtraction procedure was ap- 
plied for both reactions. A central mass band of r]^^ system near rf mass was 

2 Decay mode 77 — > 7r + 7r~7r° is not well suitable for this study because of confusion 
between charged particles originating off rf and r\ decays. 



3 



>1400 
2 1200 
^1000 
-> 800 
% 600 
^ 400 
200 





0,9 0,920,940.960,98 



m(777T + 7T ) 



GeV 



2000 
1750 
1500 
1250 
1000 
750 
500 
250 





0,9 0,9250,950,975 



m(?77T + 7T ) 



GeV 



Fig. 2. Mass spectrum of r]TT + ir~ sample produced in a) diffractive-like reaction; b) 
charge-exchange reaction. Points with errors and smooth curves represent the Real 
data and MC respectively. Vertical lines indicate the "signal mass band" for both 
samples. 

taken and also two side bands. The contribution of events from the side bands 
into distributions were subtracted from the contribution of the central band 

3 

The setup description and selection procedure for diffractive-like channel is 
described in [16]. The following selection procedure was applied for the charge- 
exchange sample. Events with two charged tracks of opposite charge and two 
showers in electromagnetic calorimeter were selected, the total event energy 
was requested at the range from 25 to 30 GeV. Events with identified e~/e + 
were excluded. Both electromagnetic showers should have energy greater than 
400 MeV, should not be associated with charged tracks. Distributions of effec- 
tive mass 771(77) are shown in Fig.l. The effective mass of the 77 pair should 
be between 455 and 620 MeV. For successful candidates, the reconstructed 
momenta of two 7 were subjected to kinematical 1-C fit to the 77 mass, and 
the fitted parameters of ?7-meson were used in further analysis. Finally, events 
with 0.9 < m(r]7[ + 7T~) < 1.0 GeV were selected. Effective mass distribution of 
(?77r + 7r — ) system is shown in Fig. 2, together with corresponding Monte-Carlo 
distribution. The fitted gaussian width of the 77' peak is ~ 7.4 MeV. 

Available statistics of 77' decays in this experiment is sufficient for precise 
determination of the Dalitz plot parameters (variations from the uniform dis- 
tribution are expected at the level of 10% ). This objective requires a careful 
Monte-Carlo (MC) simulation. The MC chain included the event generation, 

3 Mass bands of different width were tested, from 13 MeV to 21 MeV for side 
bands and the double size for the central band, respectively. The central bands of 
±21 MeV width for both reactions are indicated by lines in Fig. 2 
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Fig. 3. Distributions for the charge-exchange sample: a) correlation between mea- 
sured mass of T) meson (before fit) and its momentum in the laboratory frame for 
Real Data and Monte Carlo; b) r/ detection efficiency as a function of X; c) the 
efficiency as a function of Y variable. 



Table 1 

Dalitz Plot parameters for parametrization (3) (Diffractive-like and charge-exchange 
data separately) 



Coeff. 


Diffractive-like 


Charge-exchange 


a (at Y) 
b (at Y 2 ) 
c (at X) 
d (at X 2 ) 


-0.093 ± 0.028 ± O.OU(syst) 
-0.129 ± 0.050 ± 0.025(syst) 
-0.015 ± 0.019 ±0M8(syst) 
-0.128 ± 0.030 ± 0Ml(syst) 


-0.141 ± 0.019 ±0M0(syst) 
-0.096 ± 0.033 ± 0M6(syst) 
+0.028 ± 0.013 ± 0.020(syst) 
-0.062 ± 0.021 ± 0.009{syst) 


X 2 /ND 


69.5 / 55 


52.1 / 55 



Table 2 

Dalitz Plot parameters for parametrization (3) (combined fit) and theoretical values. 



Coeff. Combined fit Theory [13] Theory [14] Theory [15] 



a (at Y) 


-0.127 ± 0.016 ±0.008(syst) 


-0.093 


-0.116 + 0.024 


-0.132 + 0.019 


b (at Y 2 ) 


-0.106 ± 0.028 ±0.014(syst) 


-0.059 


0.000 + 0.019 


-0.108 + 0.033 


c (at X) 


+0.015 ± 0.011 ±0.014(syst) 








d (at X 2 ) 


-0.082 + 0.017 ±0.008(syst) 


-0.003 


0.016 + 0.035 


-0.046 ± 0.022 


X 2 /ND 


129.3/114 
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Fig. 4. Comparison of Real Data and Monte-Carlo weighted with fitted 
coefficients as a function of Y, in different X-intervals for charge-exchange 
reaction. 



the detector simulation based on GEANT3 package, the trigger simulation 
and the reconstruction of charged tracks. A special fast Monte-Carlo proce- 
dure was developed for gammas. Event generation for the charge-exchange 
reaction included the isotropic decay of the rj' on a background, which rises 
linearly with m(r]7m), and isotropic decay i] — > 77. For the diffract ive- like 
reaction, the distribution on mass M^^- and the polar and azimuthal angles 
of 7]' in each bin on the M v ^- were generated according to the results of Par- 
tial Wave Analysis of r)'ir~ system [17]. Reconstructed MC sample was passed 
through the same selection procedure as the real data sample. MC statistic 
contains 450000 accepted events after cuts for charge-exchange channel and 
1000000 events for diffract ive- like channel. The quality of the MC simulation 
is demonstrated in Fig.2, where the simulated mass spectra are superimposed 
on the Real Data. Another comparison is presented in Fig. 3a, where the cor- 
relation of measured 77 mass with its laboratory momentum is shown for Real 
Data and MC. Good description of this correlation in Monte-Carlo provides 
low level of systematic errors associated with 771(77) cu t. Variations of the 7/ 
detection efficiency over Dalitz plot are small, as one can judge from Fig. 3b 
and 3c. 

The fitting procedure was arranged in the following way. The Dalitz plot 
was subdivided to 8 X-bins and 8 Y-bins, i.e. 64 cells in total. The mostly 
populated cell contains ~ 250 events for the charge-exchange sample. Dalitz 
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plot parameters were obtained at the result of minimization of the functional: 

n bi „ ( r) . _ at- yr \2 

X \N U N 2 , a, b,c,d)= 2 E 1 2 j (4) 

r=l,2 i ^ri 



Here the index r enumerates two reactions considered, the index i enumerates 
cells in Dalitz plot 4 , Ni and N 2 are normalization factors, a, b, c, d are the 
Dalitz plot parameters. The M ri and D ri stand for the number of (weighted) 
entries in the i-th bin of the two-dimensional histograms in the Dalitz vari- 
ables for Monte-Carlo and for the background-subtracted data respectively. 
Statistical error a includes background subtraction and MC statistical errors. 
The MC histogram is obtained as follows: 

M ri = + aY rj + bY r ) + cX rj + dX 2 rj ), (5) 



where j enumerates the generated events and X r j and Y r j are the reconstructed 
values in the event j. 

Table 1 represents the obtained sets of parameters for both production reac- 
tions separately 5 . The quality of Monte Carlo description of Real Data is 
demonstrated in Fig. 4. Our estimation of systematic errors includes variation 
of parameters with the width of mass windows and background tuning in MC. 
Instability of the fitted parameters with respect to the variation of charge 
tracks detection efficiency in the vicinity of non-interacting beam tracks in 
the detector is small and it is also included into systematic errors. 

Table 2 contains the estimation of the Dalitz plot parameters from combined 
fit taking charge-exchange and diffractive-like channels together. The x 2 f° r 
combined fit is satisfactory, and one can conclude that 

• parameter b in parametrization (3) is negative at 3.8 stand, deviation level, 
therefore the parametrization (2) is not acceptable; 

• parameter d is negative; 

• parameter c is consistent with zero with error of order of 1%. 

Theoretical predictions from different publications [13,14,15] for parameters 
a, b, d are also given in Table 2 for comparison. 



4 Four empty cells outside the Dalitz plot boundaries are excluded. 

5 Values given for the diffractive-like channel supersede the results of previous pub- 
lication [16]. The difference is caused by more developed Monte-Carlo model in this 
study. 
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Combined fit with "linear parametrization" (2) yields unsatisfactory x 2 /ND = 
170.5 / 114 ratio and cannot be accepted. 

Attempting to clarify the origin of observed Y- and X-curvatures on the Dalitz 
plot, we adopted the procedure described in article [12] which contains a gen- 
eral Chiral Lagrangian for interactions of pseudoscalar mesons and a hypo- 
thetical nonet of light scalar mesons. However, the Y-dependence of the event 
density in this model is close to linear and cannot reproduce the observed 
quadratic term at Y parameter. 



3 Conclusions 

The Dalitz plot parameters for decay rf — > r]n + n^ in two production reactions 
are measured. We conclude that 

• results obtained for ^'-mesons produced in charge-exchange and in diffractive- 
like reactions are consistent; 

• the "linear" parametrization (2) gives unacceptable x 2 /ND which is con- 
nected with significantly negative value of the parameter b at Y 2 term in 
more general parametrization (3); 

• the value of the Charge Parity nonconservation parameter c is consistent 
with zero; 

• observed quadratic term at Y-variable is not reproduced in Effective Chiral 
Lagrangian model supplemented by nonet of scalar mesons [12], however it 
can be adopted in Chiral Unitarized model with interactions in final state 
[13,15]. The dynamical nature of this term needs clarification. 
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